Infrared and Ultraviolet QCD dynamics with quark mass for J=0,1 mesons. 
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By using a previously developed phenomenological kernel for the study of the light quark QCD 
sector and dynamical chiral symmetry breaking effects we will examine the relative infrared and 
ultraviolet QCD dynamics for J=0,1 meson properties. For the same reasons we extend and explore 
a quark mass depended generalization of the kernel in the heavy quark region and we also compare 
with the original model. The relation between the dynamics of the quark propagator and the 
effective kernel with the J=0,1 QQ and qQ mesons and quarks Compton size is also discussed. 

PACS numbers: 11.15.-q,12.38.Gc,12.38.Lg,24.85.+p 



INTRODUCTION 



A lot of effort has been focused on studying the spec- 
trum and the properties of light quark mesons (see Q- 
27]) and references therein). In these systems non- 
perturbative effects are dominant, therefore they are 
the best candidates for understanding the mechanisms 
underlying confinement and dynamical chiral symmetry 
breaking which are fundamental elements and of crucial 
importance for the theory. In some of these studies the 
quark propagator equation has provide useful insights in 
the light quark sector of QCD [§], [13] ■ It is also in- 
teresting to explore, using that fundamental block of the 
theory and the quarks bound state equation, the transi- 
tion to heavy quark physics. 

We plan to canvass, by using an effective kernel for the 
gap equation, how the quark mass affects the infrared 
and ultraviolet dressing of the propagator and how that 
in turn will alter the dynamics of the bound state of 
quarks. A more realistic case of a quark mass depen- 
dent version of the effective kernel is also explored in the 
same light. The Compton size of quarks and mesons is 
also used to qualitatively understand the relative infrared 
(IR) and ultraviolet (UV) QCD dynamics and inspired 
an approach that enabled us to reach a physical bound 
state for qb q=u/d,s,c mesons. Aspects related to the 
finite size of hadrons involving recent ideas by Brodsky 
and Shrock [ll| are also briefly discussed. 
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D^ v {k) is the renormalized dressed gluon propagator, 
T % v (p, q) is the renormalized dressed quark-gluon vertex, 
A is the regularization mass scale, with Z\, Z 2 being the 
gluon and quark propagator renormalization constants. 

Using the rainbow truncation for the gap equation and 
introducing a{q 2 ): 



Z ig 2 D^(q)Tl(p,q)^ 
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(2) 



where D l ™ e (q) is the free gluon propagator, we can disen- 
tangle the equation from the rest of the Dyson-Schwinger 
equations. 

The unrenormalized quark self-energy term of the gap 
equation in the rainbow truncation is: 



E(p)=ij{A(p*)-l} + B(p i ) = 
fA d4q G ^t m \ 

—T^kjj^Siq)^, 
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where we have set Q[k ) — Aira(k 2 ), k = p — q is the 
gluon momentum and the factor ^ comes from the trace 
over the color indexes. By taking the Dirac trace of the 
last equation we get: 



II. QUARK PROPAGATOR, MESONS BOUND 
STATE EQUATIONS AND RAINBOW-LADDER 
TRUNCATION. 



The Dyson-Schwinger equation for the quark propaga- 
tor (gap equation) has the form: 
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and if we multiply by $ and then take the Dirac trace, 
we get the second equation: 
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where we have introduced the quark propagator ampli- 
tudes (J s (q 2 ), <r v (q 2 ): 
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A{q 2 ) q 2 + M 2 {q 2 ) 
The quark propagator in terms of A , B is then: 

Z 2 (i i> + m bm ) + £ (p) = 
i ^(Z 2 + A' (p 2 ) - 1) + (m 6m + b' (p 2 )) (7) 

Using the propagator renormalization condition, 
5 _1 (p)| p 2 =Al 2 = i i/l + m r (fi 2 ), we get 

A( M 2 ,A 2 ) = 1 + A(p 2 ,A 2 ) - A' (/i 2 , A 2 ), (8) 



B{n , A 2 ) = m r ( M 2 ) + S (p , A ) - S (/r, A ). (9) 

where m r (/j 2 ) is the renormalized current quark mass 
at point /j, = 19 GeU and it is a parameter we fit to 
experimental data. 

The amplitude (BSA) Tj^(p,P) for a meson state of 
quarks of flavors a and b is given from the solution of the 
homogeneous Bethe-Salpeter equation (BSE): 
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P is the total momentum, 77 iff) is the momentum parti- 
tioning parameter for the quark (antiquark) and rj + fj = 
1, i] G [0,1]. Kt£(p,q~,P) is the unknown renormalized 
amputated irreducible quark-antiquark scattering kernel. 
Physical observables are independent of the partitioning 
parameter. 

The most general form of the BSA for psudoscalar 
mesons has four invariants while for the vector mesons 
has eight (see |U,[j|) and we use a four Chebychev poly- 
nomial expansion for each one of them. These amplitudes 
are Lorentz scalar functions of q 2 , P 2 , q.P and the mo- 
mentum partitioning parameter 77. For qQ mesons that 
parameter help us avoid having the propagator singu- 
larities inside their mass shell BSE integration domain. 
Since for the mass shell momentum P 2 — — m 2 , where m 
is the meson mass, the quark momenta in (eq. 1 10[) are 
in general complex numbers. This requires the solution 
of the gap equation in the appropriate parabolic region 
in the complex plane. 

The ladder truncation for the BSE is an approximation 
for the equation's kernel: 
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The electroweak decay constant fn of a charged pseu- 
doscalar meson [i| expressed in terms of the meson nor- 
malized BSA and quark propagators: 
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where P 2 = —m 2 H and Nc — 3 is the number of colors, 
from the trace over the color indexes. Similar expression 
exists for vector mesons. 



III. RAINBOW-LADDER EFFECTIVE 
INTERACTION. 

For the unknown effective running coupling we are go- 
ing to use a kernel-model that has been developed within 
the Rainbow- Ladder truncation of Dyson- Swchwinger 
equations. The model respects some of the most im- 
portant symmetries of QCD, like Chiral symmetry and 
Poincare covariance, while it provides quark dressing, dy- 
namical Chiral symmetry breaking, and most important, 
quark confinement. It has been used to study the physics 
of DCSB and related phenomena, like the spectrum of 
light quark mesons 

(i @, @), decay constants (0 ,0 , 
12 1, [l3j], [13]) and other physical observables (fla^ fla. 
18 1), in good agreement with experimental data ([23], 

13, m, deb- 

The so called Maris- Tandy (MT) model [4] has the 
form: 
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For the parameters we have u> = 0.4 GeV, D = 0.93 GeV 2 
and rrit = 0.5 GeV and the u/d- and s-current quark 
masses at the renormalization scale [i = 19 GeV, fitted 
to the experimental masses of pion and kaon, are m u / d = 
0.00374 GeV and m s = 0.083 GeV . For the other two 
quarks, we use the masses in fl9| . where m c = 0.88 GeV 
and m b = 3.8 GeV. The model essentially simulates 
both gluon and quark-gluon vertex dressing effects. The 
phenomenological first term defines the behavior in the 
infrared region and provides the infrared enhancement 
necessary for the right value of the quark condensate in 
the chiral limit. The second term is important in the 
ultraviolet region and is set up so that will reproduce the 
1-loop perturbative QCD running coupling behavior. 



IV. HEAVY QUARK PROPAGATOR WITHIN 
MESON DYNAMICS. 



Recently an effort has started to extend the applica- 
tions of the MT model in the heavy quarks region ([27|. 
28] 29]). Due to its behavior in the complex domain 
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FIG. 1: Re(A) for c quark in a parabolic region in the complex 
plane (P 2 = -14 GeV 2 , 77=0.50). In this plot and in the 
following ones the straight lines connecting the end-points 
have no significance and should be ignored. 



the solution of the heavy quark gap equation is challeng- 
ing. By changing the integration variable in the equa- 
tion, from the quark internal momentum to the gluon 
momentum, solved that problem and the only limitation 
comes from the propagator's singularities. The type and 
the exact location of these points is not known. Notice 
that the thranslationally invariant regularization of the 
integrals allows that change in the integration variable 
and variations of parameter rj and that our approach to 
numerically solve the equation is different than that in 
28]- 

In the case of the c quark for example with current 
mass m c = 0.88 GeV and for the parabolic region deter- 
mined by q\ — (q + r]P) 2 — q 2 - (r/M) 2 + 2ir/y / q 2 M 2 v 
where q is the BSE integration variable, with r\ — 0.50 
and P 2 = -M 2 = -14 GeV 2 (the peak of the region 
will be then at (77P) 2 = -3.5 GeV 2 ) amplitudes Re(A), 
Re(M) and Re(a s ) are plotted in Figures (jTJ [2j [3]) cor- 
respondingly. For our present studies of c-quark mesons 
we need to solve the gap equation for P 2 as small as 
P 2 = —M v ~ —9.6 GeV 2 , so the mass shell point (peak) 
is at (77P) 2 ~ —2.4 GeV 2 and the singularities are far 
from the BSE integration domain and they present no 
problem. This is also tested by varying parameter 77. 

From these plots we can see that only amplitude a s 
(and the same is true for a v ) has singularities, but not 
A and M. We can conclude then that the singularity is 
the point where the denominator of a s i v vanishes, i.e. 
q\ + M 2 {q\) = 0. 

By keeping only the infrared first term of the model 
the solution for Re(a s ) for the same quark is plotted in 
fig. ((4]). In the case where we keep only the UV term in 
the MT model, the solution has the first singularity on 
the real axis (see fig. [5]). 

Comparing the c and b quark real part of the mass 
amplitude Re(M), the last one appears to be almost flat 
near the peak of the parabolic region. Re(A) also appears 
to vary little, just above one, in the complex plane. One 
may consider that such behavior of Re(M) and Re(A) 



Re(M)(GeV) 




FIG. 2: Re(M) for c quark in a parabolic region in the com- 
plex plane (P 2 = -14 GeV 2 , 77=0. 50). 




FIG. 3: Re(a s ) for c quark in a parabolic region in the com- 
plex plane (P 2 = — 14 GeV 2 , 77=0.50) with integration over 
the gluon momentum (kcp). The peak of the parabolic re- 
gion is at (77P) 2 = -3.5 GeV 2 . We see indications of the 
existence of a pair of complex conjugate singularities near the 
peak of the region, approximately located at (x , g/ )~(-2.7, 
±3.5) GeV 2 . The type of the singularities is unknown. 
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FIG. 4: Re(a 3 ) for c quark in a parabolic region in the com- 
plex plane (P 2 = —8 GeV 2 , 77=0.50) keeping only the first 
IR term of the kernel. The peak of the parabolic region is 
at (77P) 2 = -2.0 GeV 2 . The solution is different but we still 
have the general characteristics of the solution from the com- 
plete MT model. The first pair of singularities approximately 
appear to be (x , y )~(-1.7,±2) GeV 2 . 
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FIG. 5: Re(fj s ) for the c quark in a parabolic region in the 
complex plane (P 2 — —8 GeV 2 , r/=0.50) with integration 
over the gluon momentum (kcp) and keeping only the UV 
term in the model. The peak of the parabolic region is at 
(rjP) 2 = -2.0 GeV 2 . The first singularity is on the real axis 
around (x , j/ )~(-1.9, 0) GeV 2 . 
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FIG. 6: -Re(cr s ) for the b quark in a parabolic region in the 
complex plane (P 2 — —101 GeV 2 , 77=0.50) with integration 
over the gluon momentum (kcp) with the full MT model . 
The peak of the parabolic region is at (rjP) 2 = —25.25 GeV 2 . 
Traces of a first pair of singularities can be seen near the 
peak of the region and they are approximately located at (x , 
y o )~(-240,±12) GeV 2 . 



justifies for a constituent like approximation for the b 
quark propagator. From the b quark Re(<r s ) plot we 
found the singularities have not only moved, as was ex- 
pected, deeper in the time- like region, along the real axis, 
but also further apart along the imaginary axis. This is 
exactly the opposite from a constituent-like behavior of 
a propagator. That also signifies the importance of the 
imaginary part of the dressed quark mass amplitude in 
the behavior of the propagator amplitudes a s and a v . 

In Table Q] we collect the data for the approximate lo- 
cation (x ,y ) of the singularities for the c and b prop- 
agators for the full, IR and UV gap solution. In that 
table we include a qualitative estimation of the quark 
mass dressing for each case. We assume that the loca- 
tion of the singularities along the real axis can be used 
to extract that piece of information. From these data we 



observe there is an increase in the quark mass dressing 
as we go from c to b quark using the full model. For the 
c quark there is an almost equal dressing from the IR 
and UV term of the interaction while for b quark the UV 
term provides more than twice the mass dressing of the 
IR term. 



V. RESULTS FOR QUARKONIA MESONS. 

The results for the masses and the decay constants 
of the three pseudoscalar and vector quarkonia (QQ 
Q=s,c,b) mesons are collected in Tables dTTJ) and (IIIII) . 

For the fictional pseudoscalar ss, the IR term gave a 
mass that is smaller by 18.7 % from that of the full model 
, while for the decay constant we have a 17.0 % decrease. 
For the cc meson it appears the first infrared term ac- 
counts for about 80 % of the meson mass and about 70% 
for the decay constant. Additional studies show an equal 
sensitivity of the decay constant for the ss on quark's 
UV dressing effects, directly and indirectly through the 
BSA, while for the cc it appears to be more sensitive to 
UV corrections in the meson amplitude. No bound state 
could be reached for these two systems with the UV tail 
term only. 

For the b quarkonium we first observe that now even 
the very weak, tail term, of the interaction can give us 
a bound state. Therefore a very weak attractive force is 
adequate for such heavy particles to bind them together. 
The mass of the system is lighter by just 0.055 GeV or 
0.57 % of the mass from the full model, while the decay 
constant, more sensitive to dressing details, decreases by 
12.4 %, when we keep the UV tail term only. On the other 
hand, since the IR term is much more attractive than the 
long distance term, we get a much smaller mass, a 15.4% 
relative decrease, while the decay constant decrease by 
more than 75 %. This also signifies a dominance of the 
UV tail term over the IR long distance term, supported 
by the heavy quark mass behavior of the propagator. 

Same observations are also true for the corresponding 
vector mesons data in table IIIII The only difference we 
notice is the greater gap between the full model calcu- 
lated decay constant for the b quarkonium and the one 
calculated by using only the UV term . The percent- 
age difference is about twice as much that of the corre- 
sponding pseudoscalar meson and the reason is that vec- 
tor mesons are more extended objects than pseudoscalar 
mesons, hence the IR term contributes more in the eval- 
uation of the physical observables. The decay constant 
indicates a deficiency of the low momenta behavior of the 
UV term on the quarks relative angular momentum. 

The first Chebychev moment of the dominant invari- 
ant for the equal quark pseudoscalar mesons from the 
different approaches appear in fig. ([7|). Similar behavior 
is observed for the corresponding vector mesons. 

As the mass increases we observe that the IR amplitude 
has a faster decrease with momentum than the full model 
amplitude and the system is more delocalised, consistent 



5 



TABLE I: c and b quark propagator approximate location of the singularities in the complex plane for the full MT, IR and UV 
model gap solution. All data are in GeV units. For the qualitative approximate estimation of the quark mass effective dressing 
we assume that (?7iq(19 GeV) + M®) 2 ~ —x . 
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(-1.9,0.0) 1 0.498 


b ||(-24.0,±12.0) | 


1.1 


|(-17.2,±7.5)| 0.347 1 


(-23.0,0.0)|0.996 



TABLE II: ss (fictional), cc and bb pseudoscalar meson masses 
and decay constants with their relative percentage differences 
from calculations where only the first infrared (IR) term or 
only the ultraviolet (UV) perturbative tail in the MT model 
is retained. 



1 ss,cc and bb pseudoscalar meson masses 1 


1 meson 


|| full MT||IR only|AAf/M% 


|UV only | AM/ A/% | 


| ss 


||0.696 ||0.565 | — 18.7 




1 cc 


|| 3.035 1 1 2.41 |-20.6 




\bb 


|| 9.585 || 8. 106 |-15.4 


1 9.530 |-0.6 | 


I 


cc and bb pseudoscalar meson 


decay constants | 


| meson 


1 1 full MT || IR only \Af/f% 


|UVonly|A///% | 


1 ss 


||0.182 || 0.151 |-17.0 




1 cc 


||0.387 ||0.276 | —28.7 




\bb 


|| 0.692 || 0. 172 |-75.1 


1 0.606 |-12.4 | 
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FIG. 7: First Chebychev moment of the dominant invariant 
for the ss, cc, bb pseudoscalar mesons, from the full MT model 
solution and the solution with only the IR term. For the bb 
quark system we also have an UV calculated amplitude since 
we can reach a mass shell even with only the weaker UV tail 
term in the interaction. 



TABLE III: ss, cc and bb vector meson masses and decay 
constants with their relative percentage differences from cal- 
culations where only the first infrared (IR) term or only the 
ultraviolet (UV) perturbative tail in the MT model is retained. 





ss,cc 


and bb vector meson masses I 


| meson 


|full MT| 


IR only | AM/M% \ | UV only | AM/M% \ 


| ss 


1 1.072 | 


0.949 1-11.5 || | | 


| cc 


1 3.235 | 


2.588 |-20.0 || 1 | 


\bb 


1 9.658 | 


8.130 |-15.9 || 9.586 | —0.8 | 


| ss,cc and bb vector meson decay constants | 


| meson 


|full MT| 


IRonly|A///% ||UVonly|A///% | 


| ss 


1 0.259 [ 


0.274 |+5.8 || | | 


| cc 


|0.415 | 


0.333 |-19.8 || | | 


\bb 


1 0.682 I 


0.340 |-50.1 || 0.510 |-25.2 | 



with the expectations of a relative smaller dressing of the 
quark mass provided by the IR term. The UV amplitude 
of the b quarkonium on the other hand decreases slower 
than the MT one and the system is slightly more local- 
ized. This is a sign that for this meson the UV short 



distance term is now providing most of the quark mass 
dressing and binding force for the two quarks. 



VI. RELATIVE INFRARED AND 
ULTRAVIOLET QCD DYNAMICS FOR 
QUARKONIA. 

From the last results become obvious that as the cur- 
rent quark mass increases the IR term has a marginalized 
important in the realization of meson observables. The 
question is why and how that happens. The effective 
kernel is the same and there is no change in the UV con- 
tribution in the quarks interaction. So essentially this 
is a result of the effect of the large current quark mass 
in the behavior of the propagator amplitudes and conse- 
quently the combined dynamical interplay with the MT 
effective kernel. A hrst naive qualitative explanation is 
that as we raise the quark mass the amplitudes will have 
smaller values, even very deep in the time-like region, 
and we need to get closer and closer to their singular- 
ities to notice some important increase in their values. 
That in turn will suppress the contribution and lessen 
the significance of the IR term while the UV term of the 
model will become more relevant in the evaluation of the 
observables. That is also related to the extremely fast 
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FIG. 8: Plot of the full MT model and the model with only 
the IR or UV tail term. The IR term will mostly determine 
the strong behavior of the model for k 2 < 2 GeV 2 and beyond 
that point the UV term dominates. 



decay of the IR term and much slower decrease of the 
UV one (see fig. \8§. 



A. A first look in the quark mass dynamics of the 
propagator with the MT kernel. 

To qualitatively see the effect of the heavy quark mass 
within the meson dynamics, we plot the product of the 
MT rainbow-ladder kernel with the two propagators am- 
plitudes as they appear in the BSE. There are four cases: 
a s {q\) ■ a s {q 2 _), a v (q%) ■ <J v {q 2 _) and a a {q±) ■ (?v{q%), but 
for convenience we consider only the real axis where 
we have Re(a s / v (Re(q±))) = a s / v (q 2 — (r/P) 2 ), so there 
are essentially only three amplitude products in the 
BSE integrand: a 2 s (q 2 - (ryP) 2 ), a 2 (q 2 - (nP) 2 ) and 
a s {q 2 - (nP) 2 ) ■ a v (q 2 - {nP) 2 ). The plots_of these prod- 
ucts with the MT kernel for the ss, cc, bb pscudoscalar 
mesons are in fig. dU [TUl E]) . 

We observe some mild increasing suppression in the IR 
region as the current quark mass increases for the prod- 
uct MT-a 2 , and a much stronger increasing suppression 
for the second product due to the a 2 amplitude. For the 
last product the situation is somewhere between the first 
two since it combines the effects of both propagator am- 
plitudes. For the last two cases of products involving a v , 
we actually have a small enhancement in the IR region 
for the s quark system, while from the MT-a 2 plot we 
notice a 2 has almost no effect in the IR region. For the 
UV region in all cases we have suppression of the tail 
term. That suppression though is decreasing as the cur- 
rent quark mass increases, somehow faster for the MT-a 2 
than in the MT-er s ■ a v product. Finally for the term in- 
volving a 2 there is no difference for all quark systems for 
q 2 > 10 GeV 2 (the increasing IR suppression just extends 
a little beyond 1 GeV 2 as the quark mass increase) . The 
difference in the effects between the first and the second 
type of terms are due to the fact that a s amplitude has 
also the quark mass amplitude in the numerator moder- 
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FIG. 9: Plot of the product MT-a 2 for the ss, cc, bb quark 
systems compared with the MT model behavior for external 
relative quark momentum p 2 = 0. In all cases the propagator 
amplitude along the real axis is for the on-shell solution of the 
BSE. One can clearly see the increasing suppression imposed 
by the amplitude in the IR region and the decreasing suppres- 
sion on the UV tail term, as the current quark mass increase. 
This IR suppression is mild compared to the other two cases, 
due to the presence of the dressed mass in the numerator of 
a s , moderating the decrease of the amplitude. Same things 
are true for every p 2 > but there is a significant scale down 
as p 2 increase. 
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FIG. 10: Plot of the product MT-a 2 for the ss, cc, bb quark 
systems compared with the MT model behavior for external 
momentum p 2 = . In all cases the propagator amplitude 
along the real axis is for the on-shell solution of the BSE. 
Once again we see the increasing suppression imposed by the 
amplitude in the IR region and the same, for all cases after 
certain q 2 , suppression of the UV tail term, as the current 
quark mass increase, a 2 amplitude unlike a 2 doesn't have the 
quark mass amplitude in the numerator, so we notice a much 
stronger suppression in the IR region as the current quark 
mass increase. This suppression extents somehow in the UV 
region, but for q 2 > 10 GeV 2 is the same for all cases. Finally 
notice that for the s quark system we observe an enhancing 
effect for the MT model in the IR region. Same things are 
true for every p 2 > but there is a significant scale down as 
p 2 increase. 
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FIG. 11: Plot of the product MT-a s -<j v for the ss, cc, bb quark 
systems compared with the MT model behavior for external 
momentum p 2 — 0. In all cases the propagator amplitude 
along the real axis is for the on-shell solution of the BSE. As 
one would expect, from the two previous plots, the degree of 
the effects in the IR and UV region is somewhere between the 
degree of the effects of the other two terms. Same things are 
true for every p 2 > but there is a significant scale down as 
p 2 increase. 



ating the suppressing effects. 



B. Propagator amplitudes and quark dressing 
dynamics with quark mass. 

Although from these plots one can qualitatively and 
visually understand why the IR term becomes less im- 
portant while the UV weak tail term becomes more rel- 
evant for the mesons physical observables as we increase 
the quark mass, the question is how that happens. For 
that reason we have to focus on the dynamics of the two 
amplitudes in the two regions (IR and UV) and on how 
that changes as we increase the current quark mass. The 
arguments will be again mostly qualitative. From the 
two coupled gap equations dH [5]) and using the plots 
HOP we can extract some qualitative information about 
the relative change the IR and UV tail term will bring 
upon the propagator amplitudes M(p 2 ) and A(p 2 ) as we 
increase the mass. We notice that the integrand of the 
equation for B (p 2 ) is just MT-a s . That term is almost 
unchanged in the IR region. There is only a small sup- 
pression over there as we increase the quark mass but 
there is a decreasing UV suppression. Therefore it is pos- 
sible for the UV integrated strength to dominate over the 
behavior of B (p 2 ). On the other hand MT-<7„ appears in 
the other integral equation. There is almost no change of 
that function in large momenta (larger than about 9-10 
GeV 2 ) but up to this point there is very strong suppres- 
sion greatly diminishing the contribution for the values 
of A (p 2 ). As a result we expect amplitude A(p 2 ) to 
have smaller variation in the complex plane and get val- 
ues closer to one as the quark mass increase. Amplitude 
M (M=B/A) on the other hand, since it is proportional 
to B, will still be receiving important contributions from 



the self-interaction term and even for the heavy b quark 
will vary much and not very close to the b current mass 
of 3.8 GeV. Dressing IR effects will also diminish and 
will be replaced by UV dressing effects. So it should be 
A(M E R ) < 0, A(M£ y ) > 0. If one expects that overall 
A(M E ) < then, since M s = M^ R + Mg v , we should 
have \A(Ml R )\ > \A(M% V )\. Notice that the estimated 
quark mass dressing using the singularities location of the 
IR and UV gap solutions for the c and b quarks (Table 
QJ do confirm the first two expected changes in the IR 
and UV quark mass dressing contributions as the quark 
mass is raised, but don't confirm the third case, since 
from these data appears the UV dressing effects increase 
faster than the IR decrease and the overall dressing in- 
crease. It is not clear though if this is because of a dis- 
advantage of the MT model (IR term) or because of our 
approach to estimate these dressings. From the mass am- 
plitude plots in section IV and from the quarkonia results 
in section V we have further evidence that the short dis- 
tance dressing increases in strength, faster that the long 
distance dressing decreases, with the quark mass. There- 
fore it is possible, after certain quark mass scale, that 
the heavy quarks receive more dressing from the short 
distance than the long distance self-interaction. Also, 
at a different mass scale we may also have more sort- 
distance self dressing than what the lightest u/d quarks 
get from the long distance self interaction. Another ele- 
ment that supports the validity of the above ideas is that 
the short distance analysis is essentially model indepen- 
dent, dictated by the perturbative QCD propagator and 
MT tail term. As a final note, based on the above anal- 
ysis, the heavy quark bound states and confinement of 
their quarks, are more a result of their large mass and 
their relatively larger UV dressing than a result of the 
binding interaction between them. 



C. A qualitative mathematical analysis. 

Next we will try to qualitatively understand how a 
raising quark mass can have the effects we observed so 
far. For the light quarks we assume m q << A x where 
A x ~ 1 GeV is the chiral symmetry breaking scale. . 
For the a s amplitude we have 

a s (q 2 -M 2 H /4) = 

1 M(q 2 - M|/4) 

A(q 2 - M 2 /A) [q 2 - M 2 H /A + M 2 (q 2 - M 2 H /4)] 



A(q 2 - M%/A) [q 2 - M 2 H /A + m\ + M 2 + 2m,M E ] v ' 

where Mh is the hadron mass. First we will focus on the 
IR region. We are very close to the peak of the parabolic 
region and for convenience of our analysis we take q 2 ~ 
0. For the light quarks we have strong dressing effects 
in that region so m q + Ms ~ Ms, and the propagator 
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amplitude is approximately 



My 



Ml 



Mi/4 



(15) 



Although A increases near the peak, decreasing the val- 
ues of a s , the hadron mass M# is mostly due to the 
quark IR dressing effects and therefore the difference in 
the denominator at the same time is getting smaller pro- 
viding the very small IR enhancement we noticed for the 
product of a s with the MT model. Since for the other 
amplitude the mass function is replaced by one (M — > 1) 
which is larger than m q + My; , and since the denominator 
is the same, we expect for the same reasons as before to 
have stronger than the a s IR support for the MT model. 
For the heavy quarks we assume Mh ~ 2m q . Then 



o,~i 



M| + 2m q M s ) 



(16) 



We keep the M s in the denominator because now it is 
of the order of 1 GeV. A has values close to one and 
does not affect much the behavior of the amplitude near 
the peak. Because we now have the sum of two positive 
terms in the denominator, we expect that will lightly 
decrease the values of a s near the peak. The large quark 
mass in the numerator inhibits stronger suppressing IR 
effects from this amplitude but something similar does 
not exist in <r v resulting in a stronger IR suppression. 
These observations help us to qualitatively understand 
the different IR effects we noticed in the last figures and 
how and why they change as the quark mass increase. 
We should notice at this point that the behavior of a s 
in the IR region is model dependent and for the heavy 
quarks may not present a realistic behavior, but for the 
other amplitude, a v , we know for certain there will be a 
1 jm q suppression for the heavy quarks as we increase the 
quark mass. 

The UV region analysis is easier since we are in the 
perturbative region and for both light and heavy quarks 
we have m q + My; ~ m q and A ~ 1. This is essentially 
model independent analysis. Notice at this point that 
for our analysis we are talking about dressing effects in 
the quark mass in the IR or UV region originating ei- 
ther from IR or UV interaction effects. There is though 
a shift in the momenta for the propagator amplitudes 
so in reality, when we refer to mass dressing in the IR 
region, we actually refer to the dressing in the quark 
mass near the peak (from Re(q±) = —(r/M) 2 to about 
Re(q±) ~ —(r/M) 2 + 2 GeV 2 on the real axis) of the me- 
son mass shell parabolic region in the complex plane due 
to IR and/or UV interaction effects, and when we talk 
about mass dressing in the UV region, we actually mean 
the dressing in the quark mass far from the peak (approx- 
imately for Re(q 2 S) > -(r/M) 2 + 2 GeV 2 on the real axis) 
of the meson mass shell parabolic region in the complex 
plane due again to IR and/or UV interaction effects. The 
idea of a quark receiving, at different momentum scales, 
both IR and UV dressing at the same time, can be accom- 
modated by the wave-like nature of the particle. There 



is an insignificant quark mass dressing in the UV region 
(for q 2 > 4 — 5 GeV 2 which for the propagators momenta 
in the BSE is for Re{q 2 ± ) > ~(r]M) 2 + 5 GeV 2 on the real 
axis in the corresponding parabolic region) for all quark 
masses and that's why we can assume the last approxi- 
mation for the quark mass function. The reason for this 
is that the internal quark propagator momentum in the 
gap equations depend on the external momentum p 2 and 
as the last one increase provides through the propagator 
in the gap equations integrand an increasing suppression 
on IR and UV dressing effects for the quark mass in the 
UV region. For the same reason the MT model provides 
less quark binding as p 2 increase and we move to the UV 
region, the gluon momentum in the BSE is k = p — q so 

2/2 

we have an e~ p > u term strongly suppressing IR binding 
forces in the UV region and through a more complicated 
dependence of the UV tail on p 2 we also have suppression 
of the UV binding effects again in the UV region. The 
physical interpretation is that as the quark momentum 
p 2 increase the particle will be more localized, receiving 
less IR and UV dressing and at the same time feel less of 
the color field of the other quark. We have then 



cr s (q 2 > 5 GeV 2 ) 



M H /A 



(17) 



For light quarks the numerator is very small and as q 2 
increase the amplitude will decrease very fast. For a v 
though the numerator is one which is larger than the 
quark mass and that decrease will be weaker. For the 
heavy quarks on the other hand since Mh ~ 2m q the am- 
plitude is further simplified to a s (q 2 > 5 GeV 2 ) ~ m q /q 2 . 
As the current mass further increase will provide more 
support for the values of a s in large momenta. If we 
replace that mass with one to get the other amplitude, 
since now the quark masses are larger than 1 GeV, there 
will be less UV enhancement than with the o~ s . The sim- 
ple form of the denominator also makes possible for the 
product of the two functions with the MT model to follow 
closer the changes in the behavior of the model, in fig. 
([HI ITU]) notice the obvious bending at about 2-3 GeV 2 of 
the c and b quark cases of products. Finally notice that 
we have very strong increase in the suppression effects 
in the IR region as the quark mass is raised and that af- 
fects the propagator amplitude products in a small area 
near the tip of the parabolic region, while there is a com- 
paratively mild decrease in the suppression effects in the 
UV region as the quark mass is raised in a comparatively 
much larger area of the mass shell parabolic region in 
the complex plane. Therefore, based on our analysis at 
the end of previous subsection B, the integrated strength 
of the increase in the UV dressing and binding effects 
should be slightly more that the integrated strength of 
the decrease in the IR dressing and binding effects so 
that at the end will have an overall increase in the quark 
mass UV dressing and binding energy in both IR and UV 
regions. 
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D. Physical interpretation and other aspects of the 
relative effects of the UV and IR dynamics. 

Summing up, we found that as the current quark mass 
increase the propagator amplitude a v and less o~ s will sup- 
press mass dressing and binding IR contributions while 
initially for light quarks a v will enhance IR contributions. 
At the same time o~ s will provide the support that will 
enhance UV dressing and binding effects as the current 
mass is raised. From the location of the singularities of 
the c quark propagator, from the full model calculations 
and when we keep the IR or UV term (see Table [TJ> , and 
from plots (H3 [T0| qualitatively we may conclude that this 
is the current quark mass region where this transition, in 
dominance of dressing and binding from UV over IR re- 
gion contributions, takes place. In terms of physics, the 
long distance interaction between the quarks (IR term) 
is inhibited by their large mass and now they interact 
mostly through the exchange of short-range (large mo- 
mentum) gluons. Notice at this point the unique feature 
of the theory where the interaction degrees of freedom 
(gluons) depend on the mass of their source. As the 
mass is raised becomes increasingly more difficult for the 
quarks to move further apart with the large current quark 
mass gradually replacing in that way the effect of the 
strong IR term of the interaction. In other words, the in- 
crease in the current mass now replaces the strong quark 
mass dressing IR effects of the model (enhancing at the 
same time UV dressing), making that term almost un- 
necessary, and with only a small attraction we can have 
a bound state. 

From the c quark amplitude plot in fig. Q and a sim- 
ilar one for the b quark amplitude, becomes obvious that 
the MT kernel, because of the behavior of the infrared 
term, does not support a single pole-mass constituent-like 
behavior and actually, from the b quark plot, it appears 
we are moving further away from that type of behavior 
as the current quark mass increase. We can assume then 
that a constituent-like behavior can be supported by the 
effective kernel, as we increase the quark mass, only if 
there is an important decrease in the strength of the IR 
term of the kernel. 

For the solution of the gap equation one starts by as- 
suming a free quark boundary condition at some large 
space-like momentum scale, well inside the perturbative 
region, and this solution is modified because of the in- 
teraction of the quark with the vacuum, mostly in the 
nonperturbative region for the light quarks. If we ignore 
the strong nonperturbative effects, expressed through the 
first term of the MT kernel, and keep only the pertur- 
bative UV tail term, then from perturbative QCD we 
know that the pole of the propagator, will retain some 
of its initial features determined by the current quark 
mass. The fact that we still have a first singularity on 
the real axis (fig. [5] and similar plot for the b quark) is 
a verification of that. On the other hand as we go to 
heavier quarks the IR dressing effects will become com- 
paratively smaller, the location of the singularities will 



be mostly determined by the current quark mass (and 
UV dressing) and they will be very close, along the real 
axis, to the mass-pole singularity of the free propagator. 
Therefore the dynamics of the interaction, as expressed 
through the MT model (specifically the first part) , will 
have a lesser role in the location of the singularities as the 
current quark mass is raised, but will be always responsi- 
ble for their type (confinement excludes poles) which still 
has a very subtle and important role for the calculation 
of certain meson observables. Therefore we expect that a 
small variation of the parameters of the model, will have 
insignificant impact in the location of the heavy quark 
propagator singularities. As a consequence the values 
of the physical observables will also show indifference to 
these variations and will be mostly determined by the 
quark masses. Reversing the reasoning, if the IR term 
of the model had an explicit quark mass dependence, we 
should expect a decreasing dependence of the strength of 
the model on the mass as we go to heavier quarks. 



VII. AN EFFECTIVE KERNEL WITH QUARK 
MASS DEPENDENCE. 

Higher order diagrams in the expansion of the BSE 
kernel contain internal quark propagators which will in- 
troduce a quark mass dependence to it. The parame- 
ters of an effective kernel that is applicable over a wide 
range of quark mass can be expected to have an explicit 
quark mass dependence. The way the parameters should 
vary with the quark mass is not easy to determine. In 
ref. [35[ it was noted that diagrams higher order than 
ladder-rainbow provide an attractive effect on mesons 
which decreases with increasing quark mass. Qualita- 
tive estimates of this effect for quarkonia were used to 
produce an m g -dependent effective kernel (of rainbow- 
ladder format) whose role was to produce meson bound 
state quantities that left room for subsequent and explicit 
corrections from higher order kernel processes. The IR 
strength and range of such a core kernel was fitted to re- 
produce the expected behavior of eq.(9) in ref. (3f| in the 
low quark mass region. With the MT-model form used 
as a template for the core kernel, the IR strength was 
fitted to light meson properties leaving room for higher 
order terms effects. On the other hand the IR strength 
of the MT effective kernel was fitted to light quark me- 
son properties absorbing higher order effects. Since the 
MT-model parameters are fixed these effects will be still 
present in the heavy quarks region. The extrapolation 
of the core kernel to the heavy quark region provides an 
interesting point of comparison with the MT-model. 

So far we had numerous indications that the MT ef- 
fective kernel provides too much dressing for the heavy 
quarks. At the same time we believe it gives a BSE kernel 
that is too attractive and that cancels to some degree the 
overdressing effects in the quark mass. To make things 
more complicated, we found that the heavy quark prop- 
agator itself weakens to a certain degree the model's IR 



effects. Therefore it is not so easy to determine the net 
effect of the MT interaction kernel in the heavy quark 
region from the calculated meson observables._Since the 
heavy quark region was not considered in ref. 



we ex- 



tend its application in the c and b quark region and we 
examine any benefits for studying heavy quark mesons. 
The expression for the core effective kernel of ref. (35j 

is: 



Aira{k 2 



C(lu, rh) 



(2ir) 2 k 2 



2{2n) 2 lm F{k 2 ) 



ln[r • 



fc 2 

A 5 



) 2 



(18) 



The midpoint value of u) in the minimal sensitivity re- 
gion, from studies going as far as about the s quark 
mass, was found to have the following dependence on 
the renormalization-group-invariant quark current mass 



uj(x) = 0.38 + 



0.17 

1 + x 



to 



toq = 0.12 GeV, (19) 



and in this case the product of the two core-model pa- 
rameters, D and uj, will depend only on x: 



C(x) =LuD = C a 



0.86 



1 + C 2 x + Cix 2 



(20) 



with Co = 0.11, C* 2 = 0.885, C 3 = 0.474. 

If we assume these relations are true for all current 
masses then we notice that from the chiral limit all the 
way to infinite quark mass, parameter ui(x) varies from 
w (0) = 0.55 GeV to uj(+oo) = 0.38 GeV and C(x) from 
C(0) = 0.97 GeV to C(+oo) = 0.11 GeV or for D(x) 
from D(0) = 1.763 to D(+oo) = 0.289. So the infinite 
range of quark masses is mapped into a very small domain 
of the parameters uj(x) and C(x). From the plot in fig. 
(TT2"j) we observe there is a very fast variation of C(ui) for 
quark masses up to to ~ 2 — 3 GeV and then above that 
the parameters rapidly approach their limit values and 
they don't vary much. The limiting kernel parameters 
C, u> are a possible definition of a rainbow -ladder kernel 
for very massive quarks. As a qualitative comparison 
with the MT-model at high quark mass we present some 
meson results from use of the new effective kernel which 
we call a core model. 

Since the c quark mass of 0.88 GeV at scale fi = 
19 GeV corresponds to a to of about 1 GeV, we obtain 
D = 0.308 and Q = 0.39 GeV [38] for this core kernel. 
For the b meson studies we use the core model limit val- 
ues (D=0.289, uj = 0.38 GeV). 

The results for the c- and b-quarkonia masses and their 
decay constants from the core model for comparison with 
the MT model results are collected in table (|IV[) . Sum- 
ming up, the relative percentage changes in the param- 
eters of the model and the resulting relative percentage 
changes in the quarkonia masses and decay constants, for 




0.36 0.38 0.4 0.42 0.44 0.46 0.48 0.5 0.52 0.54 0.56 

CO 

FIG. 12: Variation of parameter C( lj) = ujD{uj) vs. param- 
eter uj of the core kernel of ref. [35| for < rh < +00. The 
arrow points in the direction of increasing rh. 

the c quarkonia are: 

(AD/D,Alo/uj) = (-66.88,-2.5)% 
-> (AM/M,Af/f) Vc = (-3.3,-13.7)% 
-> (AM/M,Af/f) J/4 , = (-5.7, -23.1)% 

while for the b systems: 



(AD/D,Auj/lu) = 
+ (AM/M,Af/f) Vb 
(AM/M,Af/f) r = 



(-68.92,-5.0)% 
= (-0.50,-8.5)% 
(-0.75,-13.6)%. 



It is obvious that the decay constants, especially those 
of the vector mesons, since they are more extended par- 
ticles, are much more sensitive to parameters changes. 
Vector meson masses also appear to be relatively more 
sensitive than the pseudoscalar ones, for b quarkonia 
though the relative mass change indicate insensitivity to 
the IR dynamics of the model. Finally by comparing the 
core model results for the b quarkonium with those ob- 
tained with the UV term only in the model, we see that 
masses, especially that of the vector meson, are closer to 
the experimental ones but the vector meson decay con- 
stant at the same time becomes even smaller. 



VIII. RECOVERY OF A qQ MESON MASS 
SHELL. 



A qQ meson mass shell can not be obtained with a 
full model dressing for both quarks, and that reveals the 
limitations of the approach for light-heavy meson stud- 
ies. Given that the core model of an effective BSE kernel 
is TOg-dependent, and that its use for qQ mesons con- 
taining quarks of different masses is not defined, we de- 
cided to explore the point of view that the DSE kernel for 
the dressing of the heavy quark propagator could be too 
strong in the infrared. Hence we performed qQ meson 
calculations in which the IR term of the MT-model ker- 
nel is removed. The BSE kernel for the interaction of the 
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TABLE IV: c and b quarkonia masses and decay constants using the values from the MT and the core model [35(. The letter 
S signifies the MT-kernel values of parameters D, uj. For reasons of comparison with the core model results in the last row for 
the b quarkonia we also include the results from the UV tail term only. Experimental data are also in GeV. 



c Quarkonia 


{M™ p -,f£ p -) = (2.980,0.340) 


{Mj%,fj%) = (3.097,0.416) 


(D, cj(GeV)) 


M Vc (GeV) 


/,c(GeV) 


Mj/^GeV) 




(0.93,0.40)S 


3.035 


0.387 


3.235 


0.415 


(0.308,0.39) 


2.934 


0.334 


3.050 


0.319 


b Quarkonia 


(M^-,/^-) = (9.30,-) 


(M^ p -J^ xp -) = (9.46,0.700) 


(D, o;(GeV)) 


M Vb {GeV) 


/*(GeV) 


Mr (GeV) 


/r(GeV) 


(0.93,0.40)S 


9.585 


0.692 


9.685 


0.682 


(0.289,0.38) 


9.537 


0.633 


9.612 


0.589 


UV only 


9.530 


0.606 


9.586 


0.510 



TABLE V: qb q=u/d,s,c pseudoscalar meson masses and de- 
cay constants with their differences from experiment after 
eliminating the IR term of the MT kernel in the dressing of 
the b-quark. The full MT kernel was used for the light quark 
propagator and the solution of the BSE. All data are in GeV 
units. 



qb 


V 


M exp. 


M U H V 


AM/M% 


rexp. 
JH 


M V 
JH 


A///% 


ub 


0.90 


5.279 


4.658 


-12.16 


0.176 


0.133 


-24.4 


sb 


0.90 


5.370 


4.748 


-11.59 




0.164 




cb 


0.86 


6.286 


5.831 


-7.238 




0.453 





and mathematically accommodate the Dyson- Schwinger 
equations solution. If QCD is the correct theory for the 
interaction of quarks and gluons something like that is 
not desirable. The present studies may indicate a pos- 
sible scenario for the natural realization of these scales 
through the combined IR dynamics of the quark mass 
dependent kernel and the quark propagators in full QCD 
calculations. Some simple preliminary studies qualita- 
tively support that idea but further and more detailed 
investigation is required. 



IX. CONCLUSIONS. 



light and heavy quark though should be approximately 
unchanged since there is no significant change in the size 
of the qQ meson compared to the size of the light quark 
mesons. Solving the DSE for the b quark propagator 
with only the UV tail term of the MT kernel, where the 
full MT kernel was used for the calculation of the light 
quark propagator and the solution of the BSE, we were 
able to reach a mass shell for light-heavy mesons having 
a b quark. The masses and decay constants appear in 
Table (jV]). These calculations represent a suppression of 
IR dressing of heavy quarks in mesons. 

With this modification of the dressing of the b-quark 
we can find meson masses that are only about 7-12 % 
smaller than the experimental ones. For the only de- 
cay constant experimentally known there is a -24.4 % 
difference. The present results and the earlier results 
for the equal quark systems appear to provide a partial 
confirmation of the recent suggestion by Brodsky and 
Shrock of a maximum wavelength for quarks and gluons 

m, r 



m mesons 



11] 



7] . The existence of such a wave- 
length for quark and gluons would be due to confinement 
in mesons with size of the order of 1 fm. That in turn, 
through the Compton relation, will introduce a minimum 
quark and gluon momentum inside hadrons. Same rea- 
soning can be applied to Baryons. These scales though 
will depend on the bound state of quarks, requiring dif- 
ferent scales to be artificially introduced for each case 



The size of the QQ mesons becomes smaller as quark 
mass increases and the dressed quark quasi-particlcs 
themselves become smaller in size. The infrared compo- 
nent of the kernel relates to large distance gluons; such 
components should be less physically relevant for internal 
dynamics of heavy, small Compton size, particles. This is 
supported by our finding that the heavier QQ quarkonia 
receive diminishing contributions from the infrared com- 
ponent of our model kernel; the UV component alone 
provides a very good description of the bb states. We 
stress that this refers not only to the binding interaction 
but also to the quark self-energy dressing. 

In detail, as the quark mass increases, due to its smaller 
Compton size, UV dressing will increase and IR dressing 
will be strongly suppressed. Model independent evidence 
indicate that the short distance dressing will be, after 
certain mass scale, more than the long distance dressing, 
and at a higher mass scale, it may be over and above the 
long distance dressing of the lightest u/d quarks. The 
heavy quarks will be limited in a smaller area in space 
and the binding force will be provided by the exchange of 
sort distance gluons, signifying a smaller overall binding 
energy. The heavy quark mass and the now relatively 
stronger UV dressing, triggered by that large mass, make 
feasible, with a very weak force, to bind the particles 
together. 

On the other hand, the qQ mesons have a size that does 
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not diminish significantly with increasing heavy quark 
mass. So the infrared sector of the binding interaction 
should remain relevant. However the self-energy dressing 
of the heavy quark should not receive strong contribu- 
tions from large distance gluons. This becomes evident 
by our finding that a suppression of the infrared com- 
ponent of the b-quark dressing kernel allows a physical 
qb, q = u/d,s,c meson state even though the binding 
effective interaction and the dressing of the light quarks 
remains unchanged. 

In a similar way as the quark momentum increases the 



particle will become more localized, receiving simultane- 
ously smaller infrared and ultraviolet dressing and feeling 
less of the color field of the other quark(s) in the bound 
state. The wave-like nature of quarks can accommodate 
that interpretation. 

Aspects of this study may be related to Brodsky and 
Shrock's suggestion of a maximum wavelength of quarks 
and gluons in hadrons. The present work may actually 
suggest a possible scenario for the natural realization of 
such scales in QCD. Further studies are necessary in that 
direction. 
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